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Abstract 



CN ■ The free energy of thick center vortices is calculated in continuum Yang-Mills theory in 

> ■ 

\^ ■ one-loop approximation using the proper time regularization. The vortices are represented 



by Abelian gauge field configurations on the torus which satisfy twisted boundary conditions. 



o '. 1 Introduction 

Understanding color confinement in the framework of QCD is one of the most challenging 

p ; 

Q I non-perturbative problems of strong interaction. Center vortices offer an appealing picture of 

' confinement ^ElIHI- The vortex picture of the QCD vacuum was proposed as early as in 1978 

and has recently received support by lattice calculations. After center projection in the so 
^ ! called maximal center gauge one produces still the full string tension 0. On the other hand, 

when center vortices are removed from the Yang-Mills ensemble the string tension vanishes [H]. 
In addition, the signals of confinement disappear from QCD propagators in both the Landau 
gauge and the Coulomb gauge [Jj when the center vortices are removed. Furthermore, the 
deconfinement phase transition in the vortex picture can be understood as a depercolation 
transition (HJ and the topological charge can be understood in terms of intersections jH] and 
writhing ^H] of the center vortex sheets. 

In refs. ^2 El the free energy of center vortices has been investigated in lattice calcula- 
tions exploiting the fact that on the torus a center vortex can be induced by twisted boundary 
conditions. In ref. it was shown that in the confinement regime the free energy of a (thick) 
center vortex vanishes. In ref. ^2] the t' Hooft loop pH), which generates a center vortex [Tlj . 
has been used to calculate the free energy in the sectors of SU{2) Yang-Mills theory with 
fixed electric flux as a function of temperature and spatial volume. In the present paper we 
will carry out an analogous investigation in continuum Yang-Mills theory. We will introduce 



* supported by DFG under grant-No. DFG-Re 856/4-2 



1 



the center vortices on the torus by means of twisted boundary conditions. These boundary 
conditions are then reahzed by Abehan background fields. We calculate the free energy of 
such center vortex fields as function of the temperature in one-loop approximation^ . 

The paper is organized as follows: In the next section we summarize the relevant features of 
gauge fields on the torus satisfying twisted boundary conditions. In particular we focus on the 
derivation of the free energy of configurations with definite electric flux and its relation to the 
free energy of a static quark-antiquark pair. In section|Slwe introduce Abelian configurations 
fulfilling twisted boundary conditions. Such configurations can be interpreted as electric 
or magnetic center vortices. Then the operator of fluctuations around these gauge fields is 
determined. The spectrum of the fluctuations around such background gauge fields with non- 
zero instanton number has been calculated already in ref. pHI and in general it has negative 
modes. However, positivity of the fluctuation operator can be achieved by restriction to 
zero instanton number and by a suitable choice of some free parameters (moduli) of the 
background gauge fields. This fact, which was pointed out in ^2] for gauge fields on the 
three torus, can be trivially generalized to the four torus for gauge fields with zero Pontryagin 
index. Then the requirement of a strictly positive fluctuation operator constrains the range 
of temperatures. For the cases of purely spatial and purely temporal twists (where the 
Pontryagin index vanishes) the spectrum of the fluctuation operator is calculated. Finally, in 
sectional the proper time regularization is used to calculate the free energy of magnetic and 
electric vortices to one-loop order. Our results are then confronted with those of the lattice 
calculations of ref. |12| . 



2 QCD on the hypertorus 

To fix our notation we summarize the relevant ingredients of gauge fields on the torus T^. 
We define the four torus as R'^ modulo the lattice 

^ = |x G rV = n^L^'''^;n = (n^),n G Z^j , (1) 

where L^^^ := Lfj^e^, = 0,1,2,3 denote the vectors spanning the lattice and {e^|(e^)j/ = 
S^u} is the canonical basis of R^. Throughout this paper we choose Li = L2 = L3 = L and 
Lq = l/T with T being the temperature. Fields on are considered then as fields on R^ 
fulfilling appropriate boundary conditions. Local gauge invariants are periodic with respect 
to a shift by an arbitrary lattice vector L^^^ , whereas the gauge potential is (in general) only 
periodic up to a gauge transformation, i.e. 

Axix + l('^)) = A'l'^ix) = %{x)Axix)n^{x)^ + %ix)dxn^ix)l (2) 



^In this context let us also mention that the energy density of an infinitely thin static magnetic center vortex 
in continuum SU (2) Yang-Mills theory has been evaluated in in the Schrodinger picture to one loop order. 
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Here ^^{x) denotes the transition function in /i-direction. Transition functions in different 
directions have to respect the cocycle condition 

n^{x + = Z^Mx + , (3) 

where 

= exp -^nf,„ • Iat, n^^ G Z(mod N) (4) 



is an element of the center of the gauge group SU{N). Here we have introduced the anti- 
symmetric twist-tensors n^y and \n is the N x N unit matrix. The cocycle condition (jSJ 
expresses compatibility of two successive translations in the (/i, z^)-plane, and a non-trivial 
twist Z^i, induces a center vortex in this plane. The possible twists are divided into two 
groups. First, the spatial twists nij,{i,j = 1,2,3) which can be interpreted as the compo- 
nents of a 3-vector m = (^23, nai, 7112) which represents the direction of the magnetic flux. 
These twists induce magnetic center vortices. Second, temporal twists noi, which again can 
be interpreted as the components of a 3-vector k = (noi , no2 , Jt^os ) • These twists induce "fat" 
electric center vortices whose flux is homogeneously distributed over the whole torus. 
The Pontryagin index P]\f of a gauge potential is given by 

^iv = y d^x tr {G^^G;,) , (5) 

where 

G^, = d^A, - d,A^ + [A^, A,] = G1,Ta (6) 

is the field strength, G*^ = \€-^vai3Gai3 its dual and Ta are the generators of the Lie algebra 
of the gauge group normalized such that tr (TaT;,) = —\5ah- The Pontryagin index is fully 
determined by the transition functions 0^. In pH] it has been shown that the Pontryagin 
index is generally a fractional number 

PN = y+ (^^) ("a-) . ^ e ^ > (7) 

where v is the integer valued topological charge (instanton number) and Pf(n^i,) = ^e^uapn^unail 
is the Pfaffian of the twist tensor n^y. The fractional part of Pn is obviously due to n^y. 
Under a gauge transformation U the transition functions transform as 

^f^{x) = U{x + L^^'^)^^{x)U-\x). (8) 

The twist tensor n^j,y is invariant under gauge transformations of the transition functions. 
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2.1 Flux 



In this subsection we will derive the free energy of a pure SU{N) gauge field configuration 
on with given electric and magnetic flux (e,m). Thereby we will mainly follow the ideas 
outlined by 't Hooft [T!?] . 

The partition function under interest is given by the trace over physical states with well 
defined electric and magnetic flux (e, m). This partition function can be formally expressed 
as 

^(e, m, Lf,) := AATr {p{e, m)e-^^) , (9) 

where P(e,m) denotes the projector on definite electric and magnetic fiux, /3 = Lq is the 
inverse temperature and is a normalization constant chosen such that 3f{e = 0, m, L^) = 1. 
To obtain this projector we will first identify the physical states with given flux (e,m). To 
this end we consider the canonically quantized theory (in Weyl gauge = 0). In this gauge 
the configuration space consists of spatial gauge potentials ^j, i = 1,2,3 on the torus T'^ 
satisfying the twisted boundary conditions Q with spatial twist vector m. This vector m 
represents the magnetic fiux of the configuration under consideration. Since electric field 
operators do not commute with the gauge potential, it is harder to identify states with 
fixed electric fiux e. Physical states have to be invariant under (time-independent) gauge 
transformations 0(x). The gauge group of pure SU{N) Yang-Mills theory is SU{N)/Z{N), 
since the gauge potentials live in the adjoint representation. Therefore, gauge transformations 
take values in SU{N)/Z{N) rather than in SU{N), i.e. a gauge transformation il(x) need 
not be periodic on - but can change across the torus by a center element: 

0[k](x + L(*)) = ZiQ[k](x), i = 1,2,3, Zi:=e^\N ^ Z[N). (10) 

Correspondingly, since ■ki{SU {N) / Z {N)) = Z{N), on T'^ there are A^'^ homotopically in- 
equivalent classes of gauge transformations. The different classes are labeled by a vector 
k = {ki,k2,k^) G Z'^(mod A^) and transformations belonging to a class k are denoted by 
r2[k](x). One can choose a representative r2[k](x) from each class which takes values only in 
the Cartan subgroup of SU {N) such that 

0[ki]j7[k2] = fi[(ki + ks) mod N] . (11) 

A general gauge transformation of class k is then given by the product of r2[k](x) and a topo- 
logically trivial (i.e. truly periodic) gauge transformation (belonging to class k = 0). Physical 
states IV') have to be gauge invariant under topologically trivial (or small) gauge transforma- 
tions, but they can pick up a phase under a topologically non-trivial gauge transformation. 
Let i7[k](x) be the unitary operator which, when acting on the Hilbert space of physical 
states, generates the gauge transformation r2[k](x). Obviously, the operators 0[k](x) with 
different k commute with each other and with the Hamiltonian. The eigenstates of 0[k](x) 



4 



possess well defined electric flux as can be seen as follows. Since the r2[k](x) are unitary their 
eigenvalues are pure phase 



n[k]\iP) = e^'^('^)lV'), cT(k) SR. (12) 

and from eq. ^(T^ follows that the phase factor a is of the form 

a(k) = ^e.k (13) 

with some vector e G Z^(mod A''), i.e. we can label the eigenstate \^p) by the vector e, i.e. {ipe)- 
We will now identify the vector e as the electric flux. For this purpose we notice that the 
Wilson loop operator 

WiC) = Itr I Pexp I - ^ A^dx>^ | | (14) 

is the creation operator of electric flux lines jSDl- Let Ci, i = 1,2,3 denote a path linearly 
interpolating between points x and x + L^*). The corresponding Wilson loop operator W{Ci) 
is not invariant under homotopically non-trivial gauge transformations il[k]: 

Widf'^'''^ = ltr(0[k](x + L»)Pe-^^.^''"l7[k]t(x)) =ZiTy(Q), (15) 

where Zi is the center element defined in (^Hl- Consider now the action of 17 [k] on the state 

A[k]#(Ci)|Ve> = #(Q)^M|Ve"'^) = Z,W{Ci)nMi;,) = e^t(-k+^-')M^(Ci)|Ve) . (16) 

It is seen that the action of W{Ci) on l-i/^e) increases the i-th component of the vector e by 
one unit. Since the Wilson loop operator is the creation operator of electric flux, we have to 
identify with the electric flux vector. Accordingly the eigenvectors \'ipe) ifT^ of A[k] are 
states with definite electric flux e and the projection operator on physical states with given 
electric flux e is given by 

P(^) = ]^E^"'^'''^M- (17) 

k 

Indeed acting with P(e) on the state |V'e') one obtains: 

k k 

The desired projector P(e, m) on definite electric and magnetic flux (e, m) is then given by 
P(e) times the projector onto states with fixed spatial twist m. With the projector P(e,m) 
at our disposal the partition function for fixed electric and magnetic flux (e, m) defined by 



5 



eq. © is given by 



where 

J^k(m) = Tr (n[k]e-^^^ (20) 

is the partition function of fixed temporal twist k. (Recall that fixed spatial twist m implies 
also fixed magnetic flux.) Eq. defines the Z(A^)-Fourier transform from the temporal 
twist k to the electric flux e and is refered to as Kramers- Wannier duality. 



Formally the partition function of fixed temporal twist k (pn|l represents the thermody- 
namic average of the operator r2[k] generating gauge transformations which are periodic in 
spatial direction up to a center element, see eq. ifTT)]! . In the "coordinate" representation this 
partition function is given by (ignoring for the moment the spatial twist) 

= y" VA{A\n[k]e-^^ Pg\A) 

= J VA{A^^%~^^Pg\A) , (21) 

where Pq denotes the projector on gauge invariant states (2^. The matrix element in the 
above equation can be expressed in the standard form by a functional integral over spatial 
gauge field configurations satisfying the temporal boundary condition 

A(^,/3) = A^M(f,0). (22) 

Furthermore the projector Pq contains an integration over the gauge group with the "Haar" 
measure. This integral can be expressed as an integral over a (temporally constant) temporal 
gauge field Aq. Thereby the Haar measure becomes the Faddeev-Popov determinant in the 
gauge OqAq = 0. The resulting functional integral for the partition function is gauge 
invariant and can be expressed in an arbitrary gauge yielding 

irk= I VA^Sgf{A)e-'^^''~^\ (23) 

A(/3)=A"{0) 

where the spatial gauge field satisfies the temporally twisted boundary conditions 
Adding also spatially twisted boundary conditions introduces spatial flux m as described 
at the beginning of the section. This shows that the partition function (|2?Hl can be calculated 
from the standard functional integral supplemented by the twisted boundary conditions lfT?Hl . 
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2.2 Polyakov-loop 



The static gg-potential at finite temperature T = is defined as the alteration of the free 
energy Fqq{(3) after adding a gg-pair to the vacuum [22]. This potential can be extracted 
from the Polyakov-loop correlator 

(^(x)^t(y)^ = e-^^«^(^'y) , (24) 



where in presence of temporal twist ^q{x) the gauge invariant Polyakov loop operator is 
defined by 

x) = Itr I^Pexp J Ao(a;°,x)dxoj l^J(x)j . (25) 

Using the boundary conditions (jSJ for the gauge potential and the cocycle condition ^ one 
obtains the following periodicity properties of the Polyakov loop: 

^(x + L«) = Itr (Pe-/cfdx°Ao(xO,x+L(0)^t(^ + iW) 

= ^tr (pe-i'o'''i^''^'(^«+^o)^-ZioJ^i(x)J7j(x)17j(2; + /?)) = e~^^^!^{^) , (26) 

where ki is the i-th. component of the temporal twist vector k. After multiple usage of this 
periodicity property one arrives at 



X- Le) = e^^-''^{x) (27) 



From this relation follows that the thermic Polyakov loop correlator in presence of magnetic 
flux m is given by 

mx)^t(x-Le)) = £4?j;;§M (28) 

Comparison with eq. ipT?|l shows that ^21 

(^(x)^t(x _ Le)) = ir(e, m) = e-^^^^'"") . (29) 



Hence, on the torus the implementation of temporal twisted boundary conditions enforces 
static electric flux which can be interpreted as arising from two homogeneously but oppositely 
charged planes a distance Le apart, and thus simulates static quark and antiquark sources a 
distance Le apart (Note that {J^{'x)j3^^ {x — Le)) is independent of x). In the following we 
will realize the twisted boundary conditions by Abelian background fields. 
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3 Abelian vortex field 



In the remaining part of the paper we restrict ourselves to the gauge group SU{2). The 
implementation of flux on the torus by means of twisted boundary conditions can be most 
easily realized by Abelian background gauge fields^ of constant field strength on T^. Such a 
configuration can be interpreted as the field of a fat center vortex whose flux is homogeneously 
smeared out over the whole torus. 

The most general solution of the equations of motion of pure SU (2) Yang-Mills theory 
with constant field strength and twisted boundary conditions on reads ^H] 

Af)(x) = (^^,.x.-^)r, ^--^(q _°i)=-^'^3, (30) 

where , /i = 0, 1, 2, 3 are arbitrary constants 



is up to constant factor the field strength 

g£) = -vr^,..r (31) 

induced by the twist tensor n^^^. The gauge field (pH?)) fulfills the twisted boundary conditions 
H2|3|l with transition functions 




and twist tensor n^j^. Obviously, each non-zero component n^j^ = —Uy^ of the twist tensor 
corresponds to a non-zero field strength component G^^} of the gauge potential (EDI), which 
represents the field of \n^y\ center vortices in the (/xi/)-plane homogeneously smeared out 
on T^. We will use this field to calculate the free energy of "fat" center vortices in one- 
loop order. To this end we will compute first the spectrum and then the determinant of 
the operator of fluctuations of the gauge field around the given center vortex configuration 
lISOl- The appearance of negative eigenvalues will be avoided by appropriately choosing the 
constants Q^. 

3.1 Fluctuation operator 

The spectrum of the fluctuations around the Abelian potential IpT]! for non-zero Pontryagin 
index Q has been found already in and in general it has negative modes. On the three 
torus these negative modes can be avoided ^21 by appropriately choosing the moduli. Here we 

^Here, Abelian gauge potential means a potential that takes values only in the Cartan subalgebra of the gauge 
group. 



(moduli) and 
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will extend this result to the four torus where one has to restrict oneself to background gauge 
fields with zero Pontryagin index and to choose the constants (moduli) appropriately. 
We will shortly summarize the essential results needed below. 



The fluctuation operator is obtained by expanding the Yang-Mills action around the back- 
ground field A^*^^ IpT]) . For this purpose it is convenient to express the fluctuating gauge 
potential as 



From eqs. H2I32|) the boundary conditions for and follow: 



(33) 



6^(x + lW) =6^(x), c^(x + lW) 



exp — TTZ 



E 



c^ix). (34) 



For later use we write down the transformation properties of under a shift by an arbitrary 
lattice vector l^L^^^ which can be obtained by successive use of eq. (jSlJ: 



X + ^ lyL'^"^ = exp -m ^ Lplp^p^Xy + «vr ^ IpUpJy 

\ u / \ p,u p<u / 



(35) 



Noting that A^°^ is a solution of the equations of motion the action in terms of the fluctuations 



5Ap reads 



Sym = So 



tr -25 A,. 



6Ap{Dlf6Ap - [blSA^f + 0{6A^))] , 

(36) 



where Sq = f d^x tr | Gp'u ) is the action of the background field A^f^ and 



is the covariant derivative with respect to Ap'^ in the adjoint representation. Adopting back- 
ground gauge fixing 



A?),. 



(37) 



D°5Ap = d^6Ap+\A'^;^\5Ap 



= 



(38) 



one finds for the gauge fixed action 

S^f = So-^Jd^x {tr (SApM^^SA^) + 2tr {^Mgh^) + 0{6A^)} (39) 
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with the fluctuation operators 



M^/ = -5,,{blf - 2 [gW, •] , M,, = -{blf (40) 

and ^ being the ghost field. Using the parametrization and representing the ghost field 
as 

2 1^ 2^(4>')' -1 J 
one finally obtains up to terms of third order in the fluctuations 

(42) 

where 

Mo := {^^d^ , Mn := (^^9a - t^^XuX, + ^) . (43) 

Obviously, the operator M„ depends on the twist tensor n^^^ = —j^^- Note that the 
ghost fields (j), (/)' and r] respect the same periodicity properties as the gauge fields and 6^, 
respectively, see eq. (jHU- 

3.2 The spectrum of the fluctuation operator 

In this section we will calculate the spectrum of the operators Mq, M„ and (M„5^,y — Ani^^^), 
see eq. from which the eigenvalues of the operators Ma and Mgh l|in|l follow by using 
the decompositions (p^ll and (|^ . 

The eigenfunctions of Mq are plane waves and from the periodicity properties of the 
fields and rj one obtains the eigenvalues 

^' = Ef^)'' (44) 

These eigenvalues are also eigenvalues of Ma (|in|l (four-fold degenerate in the index /i) and 
of Mgh (non-degenerate). 

In the general case the operator Ma has negative eigenvalues implying that the considered 
background gauge potential a|{^^ is a saddle point of the action. Note that the negative 
eigenmodes occur already for covariantly constant background fields, where they signal the 
instability of the perturbative vacuum. Only if the map 

^ : ^ R^ x^^ ^^^x^ (45) 
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is degenerate (i.e. ker(.^) is non-trivial) and for suitable choice of the parameters Q^, see 
eq. JHOI), the spectrum of Af„ is strictly positive. Therefore, to ensure positivity, we will 
consider only the cases of purely spatial (magnetic) (m 7^ 0, k = 0) and purely temporal 
(m = 0, k 7^ 0) twists. In these cases the map ^ is obviously degenerate and we have the 
orthogonal decomposition^ 

= ker(^) im(^) , 
and the operator M„ decays into two parts 

each acting in one of the orthogonal spaces ker(^) and im(^). In ^n\\^cx{,9) linear term 
Tx^XyXy is absent (because (x^) € ker(^)). Therefore, the eigenfunctions of -^nlker{,_^) ^'^^ 
plane waves and can be labeled by two integers k^l. We denote the eigenvalues by \{k^i) and 
the eigenfunctions by \k^l). Eigenfunctions of can be labeled by one integer ra and 

are denoted by |m) with eigenvalue Am- Then the eigenfunctions of M„ can simply be written 
as products of eigenfunctions of -/Wnlim(^) and -^nlker(^)- 

|m, fc, I) = |m)|/c, I) =^ 
Mn\m,k,l) = \k,l) M„|ij^(^) |m) + |m) M„|^gj.(_^) \k,l) = (A^ + X{k,i)) \m,k,l) . 

The spectra of the operators Afnlim(,_^) and Afnlker(,j?) are calculated in Appendix 1X1 from 
which the spectra of the operators Ma and Mgh follow immediately and are summarized in 
tables [Hand 13 



eigenvalue 


degeneracy 


parameters 


-27r/+(p(A;,/)-q)^ 


2e 


kj ez 


2vr/ + (p(A;,0-q)' 


6e 


kj ez 


27r/(2n + l) + (p(A:,/)-q)' 


8e 


n = l,2,3,...;kj eZ 


3 2 

from eq. dH: J2 (t^) 


4 





Table 1: Spectrum of the operator M'^^ , where q :- 
are defined respectively in eqs. (f55|) . (f56|) and 



ker(,_r) 

of appendix 1X1 



. The quantities /, e and p(fc, /) 



^This decomposition is orthogonal with respect to the canonical scalar product in because the 4x4 matrix 
=^ is anti-symmetric. 
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eigenvalue 


degeneracy 


parameters 


27r/(2n + l) + (p(A;,/)- 


q) 


2e 


n = 0,l,2,...;kj eZ 


3 

from eq. (jlH): Yl 

/x=0 


I ) 


2 


1 





Table 2: Spectrum of the operator Mgh (see also the caption of table [T}. 



3.3 No twist 

For sake of completeness let us consider the case without flux m = k = 0. We will need this 
spectrum to normalize the free energy of a twist configuration with respect to the case of zero 
field strength. 

The eigenvalues of Mq are the same as in the twisted case, i.e. 



^i = Y.(-r^) ' ^e^'' ' = (^/^) (46) 

and they appear in the spectrum of M'^'^ (with degeneracy 4) and of Mgh (non-degenerate). 
The remaining spectrum consists of the eigenvalues of the operator M„ = {—id\ + 2^2a 



which read 

3 2 

Ai = X^f|^(/^-gM)') , l = {l^). (47) 

These eigenvalues appear in the spectrum of the operator M'^ (with degeneracy 8) as well 
as of the operator Mgh (with degeneracy 2). 



4 Calculation of the free energy 

In this section we will evaluate the free energy of a "fat" center vortex induced by the twisted 
boundary conditions and examine its dependence on the torus geometry, and especially on 
the temperature T = 1/Lq. For this purpose we have to evaluate the determinants of the 
fiuctuation operators. These determinants are ultraviolet singular and need regularization. 
We shall adopt the proper-time regularization which preserves gauge invariance. In the proper 
time regularization the determinant of an operator A is given by 



oo oo 



logdet^ •= ~ / — (^^P (-^^)) = - / — 5Z ^^P ' (48) 

1 ^ 1 ^ Ai5^0 



1 1 
a2" a2" 



where A is the ultraviolet cut-off. The symbol Tr' () means that the trace is taken over non- 
zero eigenvalues only. The and A^ divergencies cancel in the ratio of determinants - in 
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our case the ratio between the determinants for non-zero {n^u / 0) and zero (n^jy = 0) field 
strength. 

The proper time integral is defined only if all eigenvalues are positive. This requires 
(see the first line in table ^1 

27r/ < (p(A;,0 -q)^ VA:,ZgZ. (49) 

With given twist tensor and parameters from this inequality follows a condition for LT, 
which can be interpreted as the temperature in units of 1/L. Therefore, large values LT 
stand for high temperatures. 

After a lengthy calculation, which is performed in appendixEl one obtains the free energy 
of a "fat" center vortex in terms of the renormalized coupling constant: 



T gliT^) LT 12 LT ^2^-2 ' ^ w - , ^^m, 



+ TTT^^TTT^ log + C{fL\ Q„ LT) , (50) 



where the function C is defined in eq. (|7H|l of appendix El In the following we will separately 
discuss the cases of purely spatial (magnetic) twist m and purely temporal twist k. 



4.1 Spatial twist 

For spatial twists k = , m 7^ 0, which correspond to fat magnetic center vortices, we choose 
QtM = (IjIj^j^)- In the following table we list the values of the parameter q = — ^ 
needed in the calculation of the function C (f76)l as well as the limits for LT resulting from 
inequality (|^ : 



m 


q 


LT 


(1,0,0) 


f(Lr, 1,0,0) 


LT e IR,+ 


(1,1,0) 


f(Lr, 1,1,0) 


LT > 0.95 


(1,1,1) 


f(LT, 1,1,1) 


LT > 1.05 



(51) 



In figure^the quantity C (f76|) is plotted as a function of LT for the different values of m. One 
observes that C is nearly proportional to the number of twists (more precisely: proportional 
to ||m|p). In figure El the free energy is plotted in units of Aqcd as a function of T/Aqcd 
for the different twist configurations, where the torus extension L has been chosen to be 
L = 1/ Aqcd- The 1-loop expansion is valid if T » Aqcd = 7;, i-e. LT ^ 1. 
For large LT the free energy of thick magnetic center vortices is nearly independent of the 
temperature T. Hence, these vortices cannot be relevant for the deconfinement phase tran- 
sition. In fact it is well known that magnetic center vortices (generated here by spatial 
twists) do not contribute to the confining properties of the theory (i.e. to the temporal string 
tension). They do, however, contribute to the spatial string tension, which even slightly in- 
creases across the deconfinement phase transition. Lattice calculations [H] show that magnetic 
vortices (measured in a spatial volume at fixed time) do percolate in both the confined and 
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m = (1,(1.0) 
m = (1, 1,0) 
m = (1, 1, 1) 



10 



30 50 
LT (temperature) 



70 



Figure 1: C (fZ^ as a function of LT for different twist vectors. 



40 



0) 30 



, 20 



10 



m = (1,0,0) 
m = (1, 1.0) 
m = (1, 1, 1) 



I 1 1 1 1 1 1 

20 40 60 80 100 120 

(temperature) 

Figure 2: The free energy of thick vortices as a function of temperature. 



deconfined phase and thus cannot be used to characterize the deconfinement phase transition. 

4.2 Temporal twists 

For electric center vortices induced by temporal twists k / , m = we choose = 
(0,0, ^, \Y . In the following table the parameter q necessary for the calculation of the free 

"Here we have chosen differently from the purely spatial twist case since for k — (1, 1, 1) the vector = 
(^, ^, \, \) would be in the image of the map ^ , see eq. (QSJ, and therefore in this case we would have q = 0. 
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0.5 1 1.5 2 2.5 3 3.5 4 

LT (temperature) 



Figure 3: The quantity C (I76|l as a function of LT for different twist vectors k. The lines show the 
function C when all eigenvalues of the fluctuation spectrum are included in the corresponding proper 
time integral. The graphs plotted with the symbols (+, x, *) show C when the eigenvalues that become 
negative with rising temperature are neglected (see text). 



energy is listed. We also quote the range of temperature LT (last column) for which all 
eigenmodes of the fluctuation operator are positive, see (|i9|l: 



k 


q 


LT 


(1,0,0) 


f(0, 0,1,1) 


LT < 3.14 


(1,1,0) 


^(0,-1,1,2) 


LT < 1.67 


(1,1,1) 


i:(o,-2,i,i) 


LT < 0.60 



(52) 



In figure IHl the quantity C defined by eq. l(76)l which is part of the free energy l(5?Hl is plotted 
as a function of LT for different twists. This quantity (and consequently also the free en- 
ergy) goes to minus infinity when LT approaches the region where a single eigenvalue of the 
spectrum of the fluctuations becomes negative. This signals the limit of validity of our cal- 
culation. It is easy to see that the convergence condition can not be valid for arbitrarily 
high temperatures: The twisted boundary conditions enforce a constant flux ||k|| = JLLq 
in the plane of the twist, where flux is the product of area and (constant) field strength. 
Increasing the temperature T, i.e. decreasing the torus extension Lq, results in an increasing 
field strength until the inequality is violated. 

As can be seen from fig. 01 the free energy of electric center vortices (temporal twist) increases 
with the temperature TL. Consequently the electric center vortices become less and less 
important as the temperature increases. This result is consistent with the vanishing of the 
temporal string tension above the deconfinement phase transition observed in lattice calcu- 
lations, given the fact that electric center vortices are responsible for the temporal string 



15 



1 



0.8 - 




0.2 - 
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0.2 



(1.3 



LAqcd (torus extension) 



Figure 4: Creation probability of the thick electric vortex as function of the torus extension. 



The points (+, x,*) in figure 01 show C as a function of T for the different twists when the 
negative eigenvalues in the spectrum are neglected. The appearance of negative modes signals 
the infrared instability of the perturbative vacuum. Neglecting the negative eigenvalues the 
free energy is proportional to for large T, which corresponds to the Stefan- Boltzmann-law 
of a free Bose gas. 

Let us consider the partition function ^(m = 0) = exp { — F]^/T) of the temporal twist 
k = (1, 0, 0) on a torus with fixed ratio L/Lq = LT as a function of the temperature T. Prom 
eq. (|78|l one obtains 



i.e. for large temperature T the partition function Sf]^ decreases according to a power law 



There the partition function Sf]^ has been calculated as a function of T for different lattice 
geometries (corresponding to different values of LT) and thereby the deconfinement phase 
transition was identified. Below the critical temperature ^ is nearly one. Increasing the 
temperature ^ decreases and above the critical temperature drops to zero as predicted by 
eq. The larger LT the sharper the drop. 

The partition function in eq. can also be considered as a function of the torus 
extension L (with fixed value LT). Figure 0] shows i^k as a function^ of LAqcd for LT = 
0.3. This function can be interpreted as the creation probability of a thick vortex It 
monotonically increases with LAqcd- Although the present one-loop calculation is, strictly 
speaking, reliable only for small L it certainly shows the right tendency of for increasing 
L. Lattice calculations performed by Kovacs and Tomboulis show that the center vortex 

^In the figure L is measured in units of 1/Aqcd, i-e. with Aqcd ~ 600 MeV the value LAqcd = 0.3 corresponds 



tension. 




(53) 



with exponent —^LT. This is in qualitative agreement with the lattice results of ref. |12j . 



to L w 0.1 fm. 
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3 4 5 6 

, ^ (temperature 



, .temperature] 



(a) 



(b) 



Figure 5: Partition function of (a) temporal twist k = (1,0,0) and (b) electric flux e = (1,0,0) as 
function of tennperature. The full lines give the results under consideration of all eigenvalues. The 
plots with +, X,* show the results when the negative eigenvalues are neglected as in fig. El 



creation probability ^ = exp (— F^/T) approaches one for large L, implying that the Yang- 
Mills vacuum can be considered as a condensate of thick center vortices. 



With the above results at hand we can calculate the free energy of an electric flux e and 
the expectation value of the Polyakov-loop correlator, see eq. Since the spectrum is 

invariant under spatial rotations contributions of the twists k G {(1, 0, 0), (0, 1, 0), (0, 0, 1)} 
are identical and we will write ^k(m = 0, L^) = ^{l) if k G {(1, 0, 0), (0, 1, 0), (0, 0, 1)} and 
J^k(m = 0,L^) = ^{2) for k G {(1, 1, 0), (0, 1, 1), (1, 0, 1) and ^k(m = 0, L^) = ir(3) for 
k = (1,1,1). The partition function for the electric flux e = (1, 0, 0) is given by 

^ie,m-u,i.^j -e " 1 + 3ir(l) + 3ir(2) + ir(3) ' ^ ' 

In flgurel^the temperature dependences of iF(l) (a) and ^{e = (1,0,0)) (b) are plotted 
for different values of LAqcd- Obviously, as a function of T decreases with increasing 

temperature T and for larger L it becomes steeper at low temperatures T. In the range of 
temperature, where our calculation is valid, one observes a dual behaviour between ^(e) and 
This is expected since 3f{e) emerges from , i = 1,2,3 through a Z{2) Fourier 
transform. Qualitatively, the same dual behaviour is observed in the lattice calculation jT2]. 
Unfortunately, the perturbative ansatz has the problem that the one-loop approximation is 
valid for high temperatures. On the other hand for the calculation of 2f{e) all possible twists 
k are needed - and especially 2f{3) is known only for T/Aqcd < 0.6/{LAqcd)- Therefore, 
iF(e) can hardly be considered in the scope of the one-loop approximation. 
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5 Conclusions 



In this paper we have carried out a one loop calculation of the free energy of fat center 
vortices in SU{2) Yang-Mills theory as function of the temperature. The fat center vortices 
were induced by imposing twisted boundary conditions to the gauge field defined on a four- 
torus. These boundary conditions, in turn, were realized by Abelian background fields which 
describe fat center vortices whose flux is homogeneously distributed over the whole 4-torus. 
Accordingly the fluctuations around the Abelian background fields have to satisfy quasi pe- 
riodic boundary conditions. For arbitrary combinations of electric and magnetic center flux, 
i.e. for arbitrary twist, negative eigenmodes of the fluctuations appear. Such modes can be 
avoided for purely spatial or purely temporal twists. In the case of purely spatial twist the 
free energy is nearly independent of the temperature in agreement with lattice results |12j . 
More interesting is the case of purely temporal twist. Unfortunately in this case the range of 
validity of our calculations is rather restricted: for low temperatures the perturbative ansatz 
is not valid while for high temperatures negative modes appear in the fluctuation spectrum. 
But nevertheless in the range of validity of the present calculation our results are in qualita- 
tive agreement with the lattice results given in (T2]. Furthermore, the creation probability of 
a thick vortex in dependence on the torus extension L has been calculated and is in agreement 
with the lattice results given in 

To summarize we have been able to (at least qualitatively) reproduce the lattice results 
obtained in refs. ^21^3 by an one-loop calculation in continuum Yang-Mills theory. 
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A Spectrum of the fluctuation operators M^i and M^2 

As pointed out in section IT2l the operator see eq. (|i3|l . can be written as the sum 



Mn = M, 



+ M, 



"Iker(^) • 



In the operator M, 



"lim(^) 



the constant term 



can be absorbed in a constant shift of 



X 



2ttQx 
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The spectrum and the eigenfunctions of the operators 

2 



im(.^) 



and {M^\,^^^)6'''' -Am^n 



im(,J?) 



im(.i*' ) 



have been calculated in ^H]- In table IHl we quote the result for the spectrum, where we have 



operator 


eigenvalue 


degeneracy 


parameter 


im(.^) 


27r/(2m+ 1) 


2g 


m = 0,l,2,... 




im{,_^) 


-2nf 


2e 








2nf 


6e 








27r/(2m+ 1) 


8e 


m = 1, 2, 3, . . . 



Table 3: Spectrum of M^j^^^^) and {Mn\i^(,^)S^"' - im^^"") 



introduced the number 



/ 



|m|| ||k| 

+ 



(55) 



L2 LLo 

and the integer e which is the greatest common divisor of the integers mi, m2, ms, ki, k2, ka: 

e = gcd(mi,m2,m3,ki,k2,k3) . (56) 

Note we have chosen Li = L2 = = L and restricted ourselves to either m = or k = 0. 

The eigenfunctions of -M'nlker(^) given by plane waves exp(— ip^x^). To get its eigen- 
values one has to analyze the periodicity properties of the functions c^{x) on the space ker(^), 
or more precisely on ^\]^cv{^) ~ ^nker(^) which is a planar Bravais lattice. Let the lattice 

kcr(^) generated by the two vectors {ai,a2} (which cannot be chosen orthogonally in 
general). From eq. one obtains 



CA (x + Ziai + Z2a2) = exp j i-K^^ljuj j c\{x) , 



(57) 



where Uj = Y2^<ui^j)tJ-''^tJ-'^i^j)'^- o^*^ number the function cx is antiperiodic under 

a shift by a^, otherwise it is periodic. The possible momentum vectors can be written down 
in terms of the basis vectors {bi, b2} of the inverse lattice: 

p(A;,0 = (A; + Ai)bi + (/ + A2)b2, k,l£Z, (58) 

where Aj = if rij is even and Aj = 1/2 if nj is odd and the basis vectors hj of the inverse 
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lattice are defined by the relation 

aj • bj = 2TT5ij . (59) 

For given twist vector m or k, resp., it is straightforward to calculate the vectors slj and hj. 
The results are summarized in tabled The eigenvalues of the operator in terms 



0: 



m 


ai 


a2 


bi 


b2 


Ai 


As 


(1,0,0) 


(^0 


0,0,0) 


(0,L,0,0) 


f (1,0,0,0) 


f (0,1,0,0) 








(1,1,0) 




0,0,0) 


(0,L,L,0) 


1^(1,0,0,0) 


i(o,i,i,o) 








(1,1,1) 


{Lo 


0,0,0) 


(0,L,L,L) 


1^(1,0,0,0) 


ii(o, 1,1,1) 





1/2 



m 



k 


ai 


a2 


bi 


b2 


Ai 


A2 


(1,0,0) 


(0,0, L,0) 


(0,0, 0,L) 


f (0,0,1,0) 


f (0,0,0,1) 








(1,1,0) 


(0,0, 0,L) 


(0,L,-L,0) 


f (0,0,0,1) 


l(0, 1,-1,0) 








(1,1,1) 


(0,-L,L,0) 


(0,0, -L,L) 


|l(0,-2,l,l) 


i(o,-i,-i,2) 









Table 4: Generating vectors of ^ 



kcr(,^) 



and its dual lattice. 



of bj, Aj and Q\ read: 



\k,i) = [{k + Ai)bi + (/ + A2)b2 + 



27rQ; 



Lx 



k,l el,. 



kcr(,^)^ 



B Free energy 

To evaluate the proper time integral we start with the simplest case of zero twist n^^i, = 0. 
The corresponding eigenvalues are given in eqs. (j46|47|l . Corresponding to the different types 
of eigenvalues we will split the proper time integral log det M'^ {n^j^iy = 0) = /i + /2, where Ii 
and I2 contain all eigenvalues from equations and (|^ . respectively: 
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Under a suitable choice of all eigenvalues in I2 are positive. Using the definition of the 
theta function 



6{z, t) = exp (vrin^t + Ininz) , Im t > 0, z G C, 
and substituting r by one obtains 

00 

- logdet MAin = 0) = -Ii - Is = j ^ 



(62) 



47r 



+ 8e(r) 



(63) 



where we have introduced the function 



G(r) 



(64) 



For non-zero twist {n^^^ / 0) the calculations are done in an analogous fashion. One finds 



logdetM7(n^, /0) = 7i+/3, 



(65) 



where Ii is defined by eq. (|Hn)l and contains all the eigenvalues given in the lowest line of 
table Q, and 



h = -2 



dr 



477 



+ 2e-^^'- + 1 1 
e 7- rr^r^ 6 2 



(1 - e-fL^ 



(66) 



and 



log 



detMA(n) 
detMA(n = 0) 



00 

Is + 8 1 ^e(T) 



(67) 



47r 



In a similarly way one obtains for the determinant of Mg^: 
det Mgh{n) 



^°^det Mgh{n = 0) 
The free energy 



00 



(6S 



Fk(m) = -riogirk(m) 



(69) 
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then becomes 



i^k(m) 
T 



-So-- log det + log det Mgh 



4-7r 
l2a2 



(70) 
(71) 



where we have introduced the abbreviation 



/(r) 



■(l_e-/i^-) 
The function I(t) is regular at r = 



2e(r) 



(72) 



2\2 



/(r = 0) 



11 (/^^) 



12 LT 

so that the ultraviolet singularity in the proper time integral can be easily identified by writing 



dr 



-(/(r)-/(0)), (73) 



4tt 



4:7V 



where we introduced some small but finite Tmin > 0. The integrals on the r.h.s. are finite and 
log L^A^ is the familiar logarithmic singularity. As usual this singularity is absorbed into the 
renormalized coupling constant defined at some energy scale ^ 



1_ fx^ 
g2 + 247r2 A2 



(74) 



In the present context it is convenient to identify this scale with the temperature /i = T. 
Inserting eq. (f73]l into eq. l(7?Hl the free energy in terms of the renormalized coupling constant 
l(7^ becomes 



T 9r{T^) LT 12 LT ^ L^T^ ^ 

where we have introduced the function 



(75) 



C 



iHfL 



2\2 



12 LT 



■log- 



An 



dr 



-/(r) 



dr 



(76) 



which is independent of Tmin and can be determined numerically with very high precision. 
The function /(r) is defined in eq. ((72|l . 
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Introducing the parameter Aqcd by 



1 11 , T 



(77) 



we obtain for the free energy in terms of Aqcd- 



i^k(m,LT) 
T 



11 {fL'f 



\og{LKQCD) + C{fL\Q^,LT). 



(78) 



6 LT 
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